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We present results for the mass of the rj' meson in the continuum limit for two-flavor lattice QCD, 
calculated on the CP-PACS computer, using a renormalization gronp improved gauge action, and 
Sheikoleslami and Wohlert’s fermion action with tadpole-improved Csw Correlation fnnctions are 
measured at three values of the coupling /3 corresponding to the lattice spacing a ~ 0.22, 0.16, 0.11 fm 
and for four values of the quark mass parameter k corresponding to m-^lmp « 0.8,0.75,0.7,0.6. 

For each (/3, k) pair, 400-800 gauge configurations are used. The two-loop diagrams are evaluated 
using a noisy source method. We calculate rf propagators using both smeared and local sources, 
and find that excited state contaminations are much reduced by smearing. A full analysis for the 
smeared propagators gives m,,/ = 0.960(87)1 q’ 248 GeV in the continuum limit, where the second 
error represents the systematic uncertainty coming from varying the functional form for chiral and 
continuum extrapolations. 

PACS numbers: 11.15.Ha,12.38.Gc 


I. INTRODUCTION 

The large mass of the rf meson relative to members of 
the pseudoscalar octet has been an outstanding moblem 
in low-energy hadron spectroscopy for some time [l|,Q,|^. 
A number of lattice QCD calculations has been car¬ 
ried out II i i H i i Hi to reproduce this fea¬ 
ture, and to try to understand how it arises. These 
simulations, however, were made either with quenched 
configurations H, H, |a 0 where full QCD was em¬ 
ployed H, 0, lir)| . at a single lattice spacing. In this ar¬ 
ticle we report on a two-flavor full QCD calculation of 
the flavor singlet meson mass including the continuum 
extrapolation. The calculation is made on a set of gauge 
configurations previously generated for a study of light 
hadron physics, the results of which have been reported 
in Ref. [lH for meson and baryon spectra and in Ref. m 
for light quark masses. Three values of lattice spacing in 
the range a « 0.22 — 0.11 fm, and four values of quark 
mass covering iriTr/mp « 0.8 — 0.6 are used. 

The main computational challenge in this work lies in 
the estimation of the the double quark loop diagram con¬ 
tribution to the Tj' propagator for which the rel¬ 

ative error AGp'{t) /G{t) increases quickly with t, the 
time separation from the source. If the error becomes 
large at a time slice terr less than the first time slice of the 
plateau in the effective mass tmin, it becomes undesirable 
to fit the Gri> directly. In such circumstances the ratio 
Gp' {t) / Gtt {t) of r]' to TT propagator has been used to try 
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to cancel the effects of excited state contributions 011^. 
In this work we set out to ensure a plateau, using an 
exponential-like smeared source and the method of t/(l) 
noise source to decrease tmin- The smearing tech¬ 
nique has been previously employed in Refs. 0ll,l3,[i3- 
We also calculate the disconnected contribution with a 
local source using the method of volume source without 
gauge fixing 0. 

Preliminary results with the local source have been 
reported in Ref. 0 Here we present full results, and 
carry our analysis through in the case of the smeared 
source, where plateaux are achieved. 

The organization of this article is as follows. Details 
of numerical calculations are described in Sec. ^1 In 
Sec. Unl we give a full account of our analysis proce¬ 
dure and results including estimates of systematic errors 
arising from chiral and continuum extrapolations. Con¬ 
clusions are given in Sec. nvi 


II. NUMERICAL CALCULATIONS 
A. Action and Configuration 

We use full QCD configurations for two flavors of 
dynamical quarks. In order to reduce discretiza¬ 
tion error, these configurations are generated with a 
renormalization-group (RG) improved gauge action and 
a mean field improved clover quark action. The RG- 
improved gauge action has the form 
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TABLE I: Overview of full QCD simulations. The lattice spacing a is fixed by the vector meson mass at the physical quark 
mass and Mp = 768.4 MeV. The set of trajectories used for flavor non-singlet hadrons min and measurements of Gdisc with 
a local source m is a subset of those referred to under A^ivaj. 


p 

Csw 

LP xT 

a [fm] 

La [fm] 

K 

mps/mv 

iVlYaj 

A'skip 

-^Smeared Meas. 

1.80 

12®x24 

0.2150(22) 

0.1409 

0.807(1) 

6530 

10 

651 

1.60 

2.580(26) 

0.1430 

0.753(1) 

5240 

10 

521 




0.1445 

0.694(2) 

7350 

10 

728 




0.1464 

0.547(4) 

5250 

10 

407 

1.95 

16®x32 

0.1555(17) 

0.1375 

0.804(1) 

7000 

10 

627 

1.53 

2.489(27) 

0.1390 

0.752(1) 

7000 

10 

689 




0.1400 

0.690(1) 

7000 

10 

689 




0.1410 

0.582(3) 

5000 

10 

491 

2.10 

24®x48 

0.1076(13) 

0.1357 

0.806(1) 

4000 

5 

799 

1.47 

2.583(31) 

0.1367 

0.755(2) 

4000 

5 

776 




0.1374 

0.691(3) 

4000 

5 

767 




0.1382 

0.576(3) 

4000 

5 

785 


where are Wilson loops with size ixj, ci = —0.331, 
and Co = 1 — 8ci. For the clover quark action , we set 
the coefficient csw = where P is the plaquette 

value calculated in perturbation theory at one loop as 
P = 1 - 0.8412/3-1. 

A summary of the parameters and statistics is given in 
Table H We use three sets of configurations generated at 
bare gauge couplings (3 = 1.8, 1.95, and 2.1, correspond¬ 
ing to the lattice spacings a « 0.22, 0.16 and 0.11 fm, 
with lattice dimensions x T = 12^ x 24, 16^ x 32 and 
24^ X 48. The physical lattice sizes are roughly matched 
at La « 2.5 fm. At each /3, four hopping parameters k 
are used. They correspond to m^frup « 0.8, 0.75, 0.7 
and 0.6. The lengths of runs range from 4000 to 7350 
Hybrid Monte Carlo trajectories, and are listed under 
the column for A^Traj in Tabled 
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B. Propagator Measurements 

We calculate the single quark loop part of the ry' meson 
propagator Gconn(t) (= GT^(t)) and the two quark loop 
part Gdisc (t) both for local and smeared sources. The 
sink is always local. We use an exponential-like smearing 
kernel K 

K{\n-m\^Q) = 

a:(o) = 1 

and parameters A and B are chosen to be the same as in 
Ref. [ij • Gauge configurations are fixed to the Coulomb 
gauge. 

We try applying source smearing to both, one or nei¬ 
ther of the quark propagators. For the pion, Fig.^shows 
a case where smearing exactly one quark propagator de¬ 
livers the earliest plateau. A comparison of rj' propaga¬ 
tors for different smearings can be found in FigUl again a 
preference for smearing is generally seen, although single 
smearing is not distinguished from double smearing as it 
was in the pion case. Since this trend holds over (/3, k). 


FIG. 1: Comparison between smearing schemes for tt effective 
mass for 13 = 2.1 and k = 0.1374 on a 24® x 48 lattice. 


we focus on the single-smeared-local combination for our 
analysis. 

Our implementation of the noisy source method is to 
use C/(l) random sources, fixing to the Coulomb gauge. 
A U (1) random number exp{70(n, t)) is prepared for each 
site (n, t) of the lattice, and a smeared source is made 
according to 

ry(fi,t) = E K{n — m) (2) 

rh 

Combining the quark propagator q{n, t) for the smeared 
source r](n,t) with exp(—i6*(n, t)) yields the loop ampli¬ 
tude with a single smearing, while combining q{n, t) with 
ri'^{n,t) gives the doubly smeared amplitude. 

Our noise sources are generated only at one value of 
the (spin, color) pair at a time; other spin and color com- 
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FIG. 2: Comparison between smearing schemes for 77 ' effective 
mass for = 2.1 and k = 0.1374 on a 24^ x 48 lattice. 


FIG. 3: Relative error of one-quark smeared G'disc(t = 2) 
versus ATioise at (3 = 1.8. 


ponents of the source are left at zero. The aim of this 
procedure is to decrease fluctuations. A similar proce¬ 
dure has been developed as the “spin explicit method” 
in Ea In order to probe the gauge field evenly, we repeat 
12 times, varying the spin-color index which is chosen to 
be non-zero over all its 12 possible values, generating a 
different noisy source for each. We repeat this whole pro¬ 
cess IVnoise times, such that the total number of inversions 
for a fixed smearing is Wnv = -/Vnoise x A^spi„ x iVcoior- 

For the coarsest lattice at /3 = 1.8 the noisy source 
measurement is made with four different values of IVnoise, 
3, 5, 8 and 20. In Fig.|3we show how the error in Gdisc(i) 
varies as a function of iVnoise, taking t = 2 as a represen¬ 
tative time slice. In Fig. 01 a similar plot is shown for the 
rj' meson mass obtained by fitting the 77 ' propagator to a 
single hyperbolic cosine function. While the error is gen¬ 
erally smaller for larger Vnoise as expected, the precision 
increases only slowly with IVnoise, particularly for light 
quark masses. Furthermore these errors themselves fluc¬ 
tuate significantly. Overall, we cannot be certain that the 
Vnoise reduction will benefit our result. Since A^noise and 
computer time are linearly related, we choose Wioise = 3 
for measurements in the more costly cases of /3 = 1.95 
and 2 . 1 . At /3 = 1.8, we use iVnoise = 8 , since it gives the 
smallest observed error in the chiral limit. 

In our calculation with the local source evalu¬ 

ate G(W with the volume source method without gauge 
fixing I 3 . This measurement is made for every trajectory. 
On the other hand, Gconn is measured every Askip tra¬ 
jectories, so an additional binning of Gdisc is performed 
in preparation for constructing Gr^'. 

We calculate errors of the propagator at each time sep¬ 
aration by the jackknife method. In our study of flavor 
non-singlet hadrons El, an autocorrelation analysis has 
shown that configurations separated by 50 trajectories 
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FIG. 4: The same as Fig.0|for relative error of m^i. Data at 
the chiral limit are also shown by crosses. 


are sufficiently decorrelated, so we use bins of 50 trajec¬ 
tories in the jackknife analyses, which translates into a 
bin size of 5 (/3 = 1.8, 1.95) or 10 (/3 = 2.1) measure¬ 
ments. 

Slightly different numbers of configurations have been 
used for the rj' in the smeared case. The number of mea¬ 
surements of smeared r]' correlators for each (/3, k) pair 
is given under the column for A^smeared Meas. in Tabled 
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FIG. 5: rj' effective mass m^{t) obtained from local source FIG. 6: Ratio Gdisc(i)/G'7r(t) from local source for (3 — 2.1 
at /3 = 2.1 and k = 0.1374 on a 24® x 48 lattice. and k = 0.1374 on a 24® x 48 lattice. 


III. ANALYSIS AND RESULT 

A. Fitting of Propagators 

To extract the r]' meson mass, we fit the propagator 
using a single hyperbolic cosine function: 

Gn> (t) = Arfi [exp(-m,,/t) + exp(-m^/ (T - t))], (3) 

where T is the temporal lattice size. Alternatively one 
may fit the ratio Gdisc(t)/G,r(t) to 


= 1-B exp(-ATO t), (4) 

where Am = m^/ — m,r- The 77 ' meson mass is then 
obtained by adding to Am the pion mass m,r extracted 
from a standard fit of form 

G,r(0 = [exp (-m^t) + exp (-m^(r — t))]. (5) 

Let us first look at data obtained with the local source 
with the volume source method. Figure El shows the 
effective mass for the rj' for a typical case of (/3, k) = 
(2.1,0.1374) on a 24® x 48 lattice. The r]' effective mass 
does not show a clear plateau. We nonetheless try to fit 
the propagator for 3 < t < 6 for this example, in order to 
compare with the ratio method. Figure El shows the cor¬ 
responding ratio as a function of time. A fit of form 0 
over t > tmin yields stable values for m^' when one varies 
^min over 2 < tniin A 3. This conceals an unquantih- 
able systematic error from the effect of excited states. In 
fact, we obtain m^' = 0.670(26) from the direct ht, while 
0.584(16) from the ratio fit, showing a 14% discrepancy. 

On the other hand, fitting propagators from the 
smeared source leads to more reliable estimate of m^'. 
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FIG. 7: rj' effective mass m®?(t) from smeared source for (5 = 
2.1 and k = 0.1374 on a 24® x 48 lattice. 


We show in Fig.Qthe rj' effective mass for smeared source 
for the same simulation parameter as that for local source 
above. There is an apparent plateau starting as early as 
^min = 1- Fitting with tmin = 1 gives m,,/ = 0.528(12). 
Figure El shows the corresponding ratio Gdisc(i)/G 7 r(f) 
and the ht, which gives = 0.515(11). The difference 
of the two hts remains within 3%. 

In order to compare the quality of data from the local 
and smeared sources, we overlay two propagators and 
two effective masses in Figs. El and Eni respectively. The 
signal for t]' in the smeared case has larger errors than 
the local case. On balance, however, the advantage of 
having plateaux in the smeared case, as shown in Fig. 1101 
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FIG. 8: Ratio Gdisc(t)/G,r(t) from smeared source for fi = 2.1 
and AC = 0.1374 on a 24® x 48 lattice. 
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FIG. 9: Comparison of rj' propagator for local and smeared 
source for /3 = 2.1 and ac = 0.1374 on a 24® x 48 lattice. The 
local source propagator has been normalized to be coincident 
with the smeared source propagator at t=5. 


greatly outweighs the disadvantage of its larger statistical 
error, which can at any rate be quantified. We therefore 
concentrate on data from smeared propagators in our full 
analyses. 

Effective masses from smeared propagator at every pair 
of (/3, At) are plotted in Figs.^n(/3 = 1-8), El(/3 = 1-95) 
and El(/3 = 2.1). Consulting these effective masses, we 
determine from fitting to propagators with tmin = 1. 
Numerical values of are listed in Table QH 

For completeness, we carry out ratio fits to smeared 
propagators with tmin = 1. Numerical values are also 


^I^^I_I^I 

0 1 2 3 4 5 6 7 

t 

FIG. 10: Comparison of rj' effective mass for local and 
smeared source for jd = 2.1 and At = 0.1374 on a 24® x 48 
lattice. 
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FIG. 11: r]' effective masses m®?(t) from one-quark smeared 
sources for {3=1.% on a 12® x 24 lattice. 


given in Table HD Comparison of results from the direct 
and ratio fits are made in Fig. im We find that the dif¬ 
ference is contained within 8%. We should be aware that 
the TT effective mass does not exhibit a plateau as early as 
tmin = 1 even for the smeared source (see, e.g., Fig. 
While the amount of decrease in the pion effective mass 
is small for the smeared source, the ratio fit involves sys¬ 
tematic uncertainties of this origin. We therefore do not 
use ratio results in our final analyses. 
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TABLE II: Fitted values of m,r and for smeared source, with fit ranges. 


P K aniTT 

[^min ? ^max] 

arrir,' 

Direct Fit [tmin,tmax] Ratio Fit 

tmin 5 tmax] 

/3 = 1.8 K = 0.1409 1.155(1) 

[5,12] 

1.182(10) 

[1,4] 

1.218(8) 

[1,4] 

K = 0.1430 0.984(1) 

[6,12] 

1.019(16) 

11.4] 

1.057(11) 

[1.5] 

K = 0.1445 0.821(1) 

[6,12] 

0.921(15) 

[1.5] 

0.937(15) 

[1.5] 

K = 0.1464 0.532(2) 

[6,12] 

0.755(36) 

11,4] 

0.769(35) 

[1,5] 

b = 1.95 K = 0.1375 0.895(1) 

[7,16] 

0.960(13) 

[1,4] 

0.957(12) 

[1,4] 

K = 0.1390 0.729(1) 

[7,16] 

0.846(13) 

[1.51 

0.823(13) 

[1.5] 

K = 0.1400 0.595(1) 

[6,16] 

0.754(12) 

[1.4] 

0.716(11) 

[1.5] 

K = 0.1410 0.427(1) 

[6,16] 

0.705(24) 

11,4] 

0.653(22) 

[1,5] 

13 = 2.1 K = 0.1357 0.630(1) 

[10,24] 

0.654(9) 

[1,4] 

0.680(10) 

[1,4] 

K = 0.1367 0.516(1) 

[10,24] 

0.598(9) 

[1,4] 

0.602(9) 

[1,5] 

K = 0.1374 0.424(1) 

[10,24] 

0.528(12) 

[1.51 

0.515(11) 

[1.5] 

K = 0.1382 0.295(1) 

[10,24] 

0.450(18) 

[1,5] 

0.426(18) 

[1.5] 
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FIG. 12: r]' effective masses from one-quark smeared 

sources for /I=1.95 on a 16® x 32 lattice. 


FIG. 13: r;' effective masses from one-quark smeared 

sources for /3=2.1 on a 24® x 48 lattice. 


B. Chiral Extrapolation 

To chirally extrapolate the mass of the rj', we test 
two functional forms. One is the form appropriate to 
a Nambu-Goldstone boson, with an extra constant term 
to reflect the non-zero mass of the t]' in the chiral limit: 

{am^'Y = A{amT^Y + ( 6 ) 

We refer to this form as NGB fit. Alternatively, one may 
take T]' mass itself in a linear form in {arriT^Y, which is 
standard for vector mesons (non-NGB. fit): 

arrij^f = A{amTrY + B. (7) 

Both fitting forms reproduce our data well with compa¬ 
rable x^/d-f- = 0.9-1.1 (NGB fit) and x^/d-f- = 0.4-2.2 


(non-NGB fit), as shown in Fig's. IT?)!IT7I and ITRl We use 
the NGB fit to determine the central values, while the 
non-NGB fit is used for estimation of systematic errors. 

In order to extract at the physical point and in 
physical units, we follow the full spectrum analysis of 
Ref. and set the degenerate u and d quark masses and 
the lattice scale from the ratio of tt to p mass rriT^/mp = 
0.1757 and the p mass nip = 0.7684 GeV. 

The rj' meson mass at the physical point at each /3 
is given in Table ITTII The Non-NGB fit leads to values 
of m,,/ larger than the NGB fit. The difference is the 
largest for the coarsest lattice (16% which is about 2a) 
and smaller for the two finer lattices (8% or 1.5cr ). This 
difference yields a systematic error of about 4% (0.4cr) in 
the continuum limit, as discussed in detail in Sec. inmi 
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FIG. 14: Comparison of r/' effective mass from direct and ratio 
fit (smeared source). 



FIG. 15: Comparison of pion effective mass for local and 
smeared source for j3 = 2.1 and k = 0.1374 on a 24® x 48 
lattice. 


FIG. 16: Chiral extrapolation for /3 = 1.8 and smeared 
sources. The two extrapolated points correspond to physi¬ 
cal and to zero pion mass. 
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FIG. 17: Chiral extrapolation for /3 = 1.95 and smeared 
sources. The meaning of symbols is the same as in Fig. Ei 


C. Continuum Extrapolation 

FigtirelT^shows mr^i as a function of a. We use a linear 
form for the continuum extrapolation, 

rrijji = C + Da (8) 

to estimate the central value in the continuum limit, since 
we employ a tadpole-improved value for the clover coef¬ 
ficient in our quark action. The y^/d./. of the fit is 4.2. 
We also try a constant plus quadratic form, since one 
may expect 0(a) effects which remain after tadpole im¬ 
provement are small: 

rrin' = C + Da^. (9) 


We find x^/d./. of this fit to be 2.8. Finally, since the 
data hardly change between the finest two lattice spac- 
ings, we try removing the coarsest point and fitting to 
a constant. The quadratic and constant fits are used to 
estimate the systematic error. 


D. Systematic Error Estimate and Final Result 

We now consider the systematic errors based on the 
two variant forms for the continuum extrapolation, and 
the alternative form for the chiral extrapolation. In Ta- 
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FIG. 18: Chiral extrapolation for /3 = 2.1, smeared sources. FIG. 19: Continuum extrapolation of 77 ' meson mass obtained 
The meaning of symbols is the same as in Fig. uni with smeared sources. 


TABLE III: Comparison of 77 ' mass at the physical point from 
Eqs. 0 (NGB fit) and lO (non-NGB fit). 


P m„,{GeV) 

NGB non-NGB 

1.8 0.509(39) 0.589(24) 
1.95 0.714(36) 0.766(26) 
2.1 0.709(41) 0.764(30) 


bleCVI we list the four estimates of masses in the con¬ 
tinuum limit. The central value is obtained from the 
linear continuum extrapolation of results from the NGB 
fit. The quadratic continuum extrapolation gives a lower 
estimate, and the constant extrapolation a still lower one. 
On the other hand, the non-NGB chiral fit gives a raised 
estimate. We therefore take the difference between the 
central value and the value from the constant continuum 
extrapolation as the lower part of systematic error, and 
the deviation of the non-NGB chiral fit as the upper part. 

Our final result for rriri' in the continuum limit reads 

run' = 0.960(87)ti^;^i^^ GeV. (10) 

IV. CONCLUSIONS 

We have made the first calculation of the mass for the 
flavor singlet pseudoscalar meson in the continuum limit. 
We find that, using smearing, it becomes possible to fit 
the flavor singlet pseudoscalar meson correlator directly 
with a hyperbolic cosine ansatz. Similar conclusions have 
been reached in Refs. [1113 

There is a lower systematic uncertainty of nearly 30% 
in our final result. The systematic error breakdown indi¬ 


TABLE IV: Systematic variations in tti,,/ over fit forms for 
chiral and continunm extrapolations. 


_m„/ (GeV) 

central value 0.960(87) 

quadratic continuum extrapolation 0.819(50) 
constant continuum extrapolation 0.712(27) 
non-N. G. B. chiral extrapolation 0.997(61) 


cates control over the continuum extrapolation to be the 
most important aim for future simulations, although this 
control may be established indirectly by some pattern of 
reduction in statistical error. On the other hand, the 
higher systematic error, coming entirely from the chiral 
extrapolation, is only 4%. 

Our continuum result is in agreement with the experi¬ 
mental value for the rj' mass of rur,' = 0.956 GeV, despite 
the fact that our calculation is carried out for a flavor sin¬ 
glet meson composed of u and d quarks and within two- 
flavor approxaimation to QCD. Thus, calculations are 
underway to attempt to extract the elements of the mix¬ 
ing matrix between quark-based states and eigenstates 
of mass in the framework of the three-flavor QCD with a 
partially quenched strange quark. We also aim to follow 
this analysis with a continuum result from “2-1-1” flavor 
QCD, in which the strange quark is also treated dynam¬ 
ically, with similar or better statistics than the present 
work. 
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